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I. INTRODUCTION 

One of unexpected experimental discoveries in heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) is 
that the created matter could be close to a perfect fiuid[l[. Then, considerable changes have occurred with our under- 
standing of quark-gluon plasma (QGP), which is now believed to be rather strongly coupled or strongly interacting 
just above the critical temperature (T c ) of the chiral and deconfmcmcnt transitions, T > T c . Nonperturbative results 
from lattice QCD seem to support this view: The lowest charmonium state survives for T > T c , which suggests strong 
correlations between quarks [2|. Furthermore, very small viscosity-entropy ratio suggested from RHIC experiments is 
consistent with the result of the strongly coupled supersymmetric Yang-Mills theory based on the analysis with the 
AdS/CFT correspondence [3|. 

The present paper is concerned with spectral properties of fermions in such a strongly coupled system, because 
quarks are the basic degrees freedom of the created matter in the deconfined phase together with gluons at RHIC. 
It is not even clear whether quarks manifest themselves as well-defined quasiparticlcs in such a system. So far, there 
exist some theoretical investigations on this issued, Q. Recently, the thermal mass of the quark near T c was studied 
in quenched lattice QCD and was fitted well with a two-pole ansatzQ. The authors and Yoshimoto investigated the 
quasiquark spectrum with a strong coupling gauge theory based on the Schwinger-Dyson equation (SDE) and found 
that the thermal mass is of the order T and depends on the coupling little in the strong coupling region 0. In Ref. [|, 
it is shown in the Nambu-Jona-Lasinio model that the quark spectrum has three peaks near T c , which result from the 
coupling with fluctuations of the chiral condensate Q. 

The above peculiar three-peak structure of the quark spectrum is then found to be rather generic at finite T: It 
is due to the coupling with massive bosonic modes through the Landau damping, which was analyzed in detail in 
Yukawa models at one loopfiol Ws\. While the three-peak structure appears most clearly for the massless fermion, it 
even appears for small massive fermions [ill]. One of the remained issues is whether it is stable against higher order 
effects. One of the purposes of this paper is to study it in a linear sigma model by incorporating nonperturbative 
effects with the SDE[l2j]. It is shown in Ref. [l(| that for the massless fermion, the three-peak structure appears 
independently of the coupling over some range of T at the one-loop order. Because the SDE should reproduce the 
perturbation theory in the weak coupling region, the three-peak structure would remain there even in this approach. 
However, it is quite nontrivial whether it still appears in the strong coupling region where higher order effects are not 
negligible. 

Another purpose of this study is to formulate coupled SDE for fermion and boson at finite T. In t he p ast approaches, 
the boson propagator was not solved self-consistently, but fixed at the beginning 
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far as we know, this paper gives the first numerical result for the coupled SDE for a fermion-boson system at finite 
T. The point in our formulation is that while the imaginary parts of the self-energies are evaluated through loop 
integrals, the real parts are evaluated from the imaginary parts with the dispersion relation which is a one-dimensional 
integral, and thus the computation is much faster than that of the four-dimensional loop integral. This reduction of 
the computational cost makes it possible to solve the coupled SDE for fermion and boson quite fast. 

The content of this paper is as follows. We employ a linear sigma model in the chiral symmetric phase to study 
the fermion spectrum at finite T . In Sec. [Til the model is introduced and the SDE arc formulated using the real 
time formalism. We use the ladder approximation to construct the SDE in which the bare vertex is used for the self- 
energies. Then, the renormalization is carried out for the boson mass and some remarks are made for the triviality 
of the model. In Sec. IIIH some numerical results of the fermion spectral functions are shown, focusing on the T 
and coupling dependences and nonperturbative effects by comparing with the results in the one-loop calculation. 
Concluding remarks and outlook are given in Sec. IIVI 
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FIG. 1: Self-energies for fermion (left) and boson (right). Solid lines denote the fermion propagator and dashed lines the boson 
propagator. 

II. SCHWINGER-DYSON EQUATIONS AT FINITE TEMPERATURE 

We consider a U(l)iX U(l)^ invariant linear sigma model with Dirac fermions, 

C = + -^ikP^M - ^M*M + g$ L M$ R + $ R M*il> L ], (2.1) 

where M = <f>i +i<f>2 and (j>\{_4>2) denotes a scalar (pseudoscalar) boson field. Because our motivation lies in the quark 
spectrum in the QGP phase with possible bosonic excitations, we assume m§ > 0, i.e. the chiral symmetry is not 
spontaneously broken. Self-coupling terms of the boson fields are not introduced for simplicity. Their main effect is to 
shift the thermal mass and the width of the boson spectrum. As mentioned below, we renormalize the thermal mass 
of the boson and measure all the dimension-full quantities in units of the renormalized boson mass. Then, effects of 
the self-coupling of the bosons are not important for the fermion spectrum. Although we consider a chiral U(l)i x 
U(1)_r symmetry for simplicity, an extension to the SU(2)lX SU(2)a symmetry is trivial as mentioned below also. 

A. Formulation of the Schwinger-Dyson equations 

To construct the SDE, we employ the real time formalism, ft is convenient to use the following Green functions, 

iS > (x,y) = (ip(x)ip(y)), iS < (x,y) = (-ip(y)ip(x)}, 

iG>{x,y) = {<j )i {x)4> i [ 3 j)), iGf{x,y) = {4> i {y)4> i {x)), (i = 1,2). (2.2) 

Because we consider the chiral symmetric phase, there is no difference between Gf and GJ, and thus we simply write 
them as in the following. 

We formulate the SDE in the ladder approximation whose self-energies are depicted in Fig. [TJ The fermion self- 
energy reads 

= e c (x° - y°)S > (x, 2/ ) - 9 c (if - x°)X<(x, y), (2.3) 

where 

T&(x,y) = -2g 2 iSZ(x,y)G*(x,y), (2.4) 

and 9 C denotes the step function along the closed time path in the real time formalism[2l|. The factor 2 in Eq. (|2.4p 
comes from the number of boson fields, 4>\ and (j> 2 . Thus, for the SU(2)^x SU(2)# invariant linear sigma model, this 
factor is replaced by 4, which is the only difference between the U(1)l X U(l)# model and the SU(2)l X SU(2)^ model 
in the fermion self-energy. The boson self-energy reads 

U(x,y) =9 c (x°-y°)Il>(x,y)+8 c (y -x )n<(x : y), (2.5) 

where 

n*(a;,y) = -g 2 Wr[S*(x,y)S*(y,x)], (2.6) 

and Tr means the trace for the Dirac index. For the SU(2)lX SU(2)a invariant model, a factor 2 coming from the 
trace in the flavor index is multiplied in Eq. (|2.6|) . 
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In equilibrium, we can Fourier-transform the above functions using the translation invariance. In momentum 
space, the Green functions S<(G<) have simple expressions with the spectral function Pf(pb) and the Fermi-Dirac 
(Bose-Einstein) distribution function f(n), 

iS>(p°,p) = [1 - /(po)]p/(po,p), iS<( P °,p) = f(p )p f (p°,p), (2.7) 
iG>(p ,p) = [l+n(j> )]pb(p ,p), iG < (p ,p) = n(p )p b (p°,p). (2.8) 
Using these expressions, the self-energies in momentum space are written as 
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(2.10) 

with P = (p°,p), Q = (<?°,q), Pj = (p°,Pf) and pf — 70/9° — 7 • Pj. The functions £ p and n p are related to the 
corresponding retarded self-energies, 

Im£*(P) = - is p (P), Re£«(P) = P /" ^^^Vp) 



ImII*(P) = -in p (P), ReIT*(P) = P f (2.11) 



in„(P,, R en«(P, = p/^^ 

where P denotes the principal integral. 

The retarded Green functions are expressed with the retarded self-energies and thus S p and II p as 

S R (P) = ' 



If + ReS«(P) + ilmE R {P) ' 

G ^ (P) = P 2 - mg + Rcn«(P) + Tmn fl (P) ' (2 ' 12) 

The spectral functions are in turn expressed with the retarded Green functions, 

Pf (P) = -2lmS R (P), Pb {P) = -2ImG J? (P). (2.13) 

Equations (|2.9p - (|2.13p are closed form for the spectral functions pf and p b , and thus give the self-consistency condition. 
Two comments are in order. 

One is that this formulation is similar to that given in Appendix D of Ref.[22j. The main difference is in the 
calculation method of the self-energies: While we evaluate them as momentum integrals, in Ref.[13| the self-energies 
are Fourier-transformed and evaluated in the coordinate space. In the latter, the fast Fourier transformation algorithm 
can be used for the numerical calculation and in general is faster than the direct evaluation of the four-dimensional 
integral. In equilibrium, however, one of the angle integrals is trivial and thus the direct evaluation of the integrals is 
effective as well. 

The other is that in this method, only the imaginary parts of the self-energies are evaluated through S p and II p . 
The real parts are obtained from the imaginary parts by the Kramers-Kronig relation which is given by the one- 
dimensional integral. Therefore, this is much efficient than the method in which both the real and imaginary parts of 
the retarded self-energies are evaluated with loop integrals. This is the main reason why we can compute the coupled 
SDE at finite T with a less numerical cost. 



B. Renormalization 



In Fig. [I] the boson self-energy is quadratically divergent and the fermion self-energy is logarithmically divergent, 
and thus both have to be regularized. In the numerical calculation, we introduce an ultraviolet cutoff in the momentum 
integrals. While the linear sigma model is renormalizable perturbatively, it is known that there is no nontrivial con- 
tinuum limit. Here, we do not pursue this issue and the renormalization is carried out partially: We only renormalize 
the quadratically divergent boson self-energy using the once subtracted dispersion relation, 



Rett R < 



{Po,P )=c+ (Po-a )P / — - 2 £-r^ 2V 2 ' 14 ) 

Jo 2tt (pi - q£)(q£ - a 2 ) 



4 



where p = |p| and we have used a property that Tl p (po,p) is an odd function for po. The subtraction constant c and 
the renormalization point a are determined by imposing the 'on-shell' renormalization condition, 



c = 0, a = \J 



P 



mo, 



(2.15) 



where 'on-shell' means that the renormalization is done so that the boson has the dispersion relation po = ^Jp 2 + rn 2 
including thermal effects. In other words, we renormalize the boson mass such that the thermal mass of the boson 
should be m. This T-dependent renormalization is contrasted to the T-independent renormalization used in the one- 
loop calculation where the T-independent terms can be separated from the others and are only renormalized [Icl [TlT| . 
In the Schwinger-Dyson (SD) approach, however, the T-independent and T-dependent terms cannot be separated 
from each other by the self-consistency condition and thus such a renormalization is hard to impose. Because we are 
mainly interested in the fermion spectrum, we employ this simple 'on-shell' renormalization for the boson mass. 

The boson self-energy is still logarithmically divergent after the above renormalization and the fermion self-energy 
is also logarithmically divergent, both of which are in principle renormalized with the wave function renormalization. 
We, however, leave an ultraviolet three-momentum cutoff as in the previous workQ. One reason is that numerical 
calculations do not suffer from such a logarithmic divergence very much. The other is that the wave function renor- 
malization modifies only the over-all factor in the propagator, which may not change the relative peak structure of 
the spectral function. Our main interest in this paper is to investigate the fermion spectrum at low energy and low 
momentum where thermal effects rather than quantum effects are dominant. 

Therefore, the self-energies are written with the cutoff as follows: 
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2tt 3 

9' 
4ir 3 p 



g 2 f A f A 

s °(i?o,0) = -^3 / dq / dqq 2 [l- f(q Q )+n{p ~qo)}po(qo,q)Pb(po-qo,q), 
-A Jo 

,2 fA /-A pp+q 

dq dq dkqk[l- f(q ) +n(p - q )}p (q ,q)pb(po- qo,k), 

-A Jo J\p-q\ 



E p y(p ,0) = 0, 
Z P v(po,p) -- 



B7T 3 p 2 J_ A 
n 2 r A 



dqa / dq 



p+q 



n P (p ,p) = - 



2ir 3 p 



A rA rP+q 
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qkp (q ,q)po(po - qo,k) Pv(qo, q)pv(Po ~ 9o, k) 



(2.16) 

(2.17) 
(2.18) 



dkk(p 2 + q 2 - k 2 ) [1 - f(q ) + n(p - q )] pv(qo, q)Pb(po ~ qo, k), (2.19) 



n p (po,0) = -Hr / dq dqq [l-f{q )-f{po-qo)][po(qa,q)po{po-qo,q)+Pv(qo,q)pv(po-qo,q)], (2-20) 



(2.21) 



with T, p {po,p) = 7 Sp(p ,p) - 7 ' P^ p v(po,p), Pf{Po,p) = l°Po{Po,p) ~ 1 ■ PPv(po,p) and p = p/p. It is noted 
that more strictly, the cutoff A is introduced so that the self-energies T, p (po,p), U p (po,p) and the spectral functions 
Po,v,b(po,p) have nonzero values in the range —A < po < A and < p < A. 

As mentioned above, while the boson mass is renormalized, the coupling constant g is not and thus depends on the 
cutoff A. In the numerical calculation given in the next section, we take A/mo ~ 10 which is determined from the fact 
that larger values of A are hard to calculate for the numerical accuracy and smaller values of A would affect thermal 
effects significantly because we set T/mo ^ 2. 1 We have checked that the following results of the fermion spectrum 
are stable under small change of the values of A, because we focus on the energy-momentum region where thermal 
effects such as the Landau damping are dominant. 



As mentioned in the main text, the triviality means the coupling g at infrared goes to zero for A — * oo. However, if we focus on the 
spectrum for po,p,m T, the temperature would be an effective ultraviolet cutoff, and thus the system would be renormalizable even 
nonperturbatively owing to effective three-dimensional reduction. Although we consider the spectrum up to T ps m, we hope that such 
an effective reduction partly works. 
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FIG. 2: The fermion spectral functions p+m for g 2 = 1 and T/m = 0.6 (top left), 0.8 (top middle), 1.0 (top right), 1.2 (bottom 
left), 1.6 (bottom middle), 2.0 (bottom right). The figure is clipped at p+m = 50. 



III. NUMERICAL RESULTS 



In this section, we examine how the quasiparticle picture of the fermion changes by studying the fermion spectral 
function as temperature or the coupling changes. It is useful in the following analysis to employ the expression 



P±(po,P) = Po(Po,p) ± Pv(po,p), 



(3.1) 



where p + and p~ represent the spectrum of the fermion and the antifermion excitations, respectively, and are obtained 
from the fermion self-energies [ToJ as 



P±(po,p) = -21m 



PoTP- S±(p ,p) 



(3.2) 



with £± = Sg- =f Ey. Because we have a relation p_(po,p) = P+{~ POiP) by the charge conjugation invariance, we 
show only the numerical results for p + below. 



A. Temperature dependence 

Before showing the numerical results, we briefly review the fermion spectrum at T — and high-T limit in our 
model. At T = 0, the renormalization prescription for the boson mass given in Sec. Ill Bl becomes exactly the on-shell 
renormalization, and thus the poles of the fermion propagator are on the light cone. The peaks of the fermion spectral 
functions p + and p_ are then along the light cone, p = p and p = —p, respectively. Effects of the interaction through 
the SDE arc contained in the widths of the peaks and the continuum part. 

In the high-T, weak coupling limit, the fermion spectrum approaches that calculated in the HTL approximation 
with the thermal mass of the fermion tut — v / 2oT'/ 4[23j . Both p + and p- have two 5-functions corresponding to 
the normal quasiparticle and plasmino excitations, in addition to the continuum part in the spacelike region. In the 
high-T limit but with fixed coupling g, these two peaks have widths of the order of arriT [2l[. Even in the SD approach, 
the fermion spectrum should approach that in the HTL approximation for the weak coupling because the self-energies 
shown in Fig. Q] include all the hard thermal loops of this model and the perturbation theory should work. 

In Fig. [21 we plot the fermion spectral function p + for various temperatures. At T/m = 0.6, the peak position of 
the spectral function p + is almost on the light cone; p+ has sharp quasiparticle peaks around pa = p. At T/m = 0.8, 
we see that the quasiparticle peaks are deformed near p/m — 0.2 and that the peak position deviates from the light 
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FIG. 3: Temperature-dependence of the fermion spectral function p+m for p = 0, g = 1. 
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FIG. 4: The fermion self-energy E+/m for p = 0, g 2 = 1: the imaginary part (left panel) and the real part (right panel). 



cone for p/m < 0.2. This deformation is more remarkable as T is higher: At T/ra — 1.0 and 1.2, the quasiparticle 
peaks start to split for p/m < 0.2, and there appear broad peaks both in the positive and negative energy regions 
around p = 0. Therefore, we see a three-peak structure in the low- momentum region. At T/rn = 1.6, three peaks 
clearly separate from one another and the strength of the central peak becomes weaker than those of the peaks in 
both sides. At T/m = 2, the central peak almost disappears and thus only two peaks remain in the low-momentum 
region. These two peaks correspond to the normal quasifermion and the antiplasmino seen in the spectrum in the 
HTL approximation (2^| . The fact that these peaks survive for the strong coupling region is consistent with that in 
Refs.0,121- 

In the same figure, as momentum is higher, the peak of the normal quasifermion becomes larger and sharper. This 
is because thermal effects get weaker and thus the spectrum should approach the free one up to the width due to 
quantum corrections as momentum is higher. The antiplasmino peak in turn rapidly suppresses as momentum is 
higher, which is consistent with the behavior of the HTL approximation where the residue of the (anti-)plasmino peak 
suppresses exponentially with momentum. 

The above result resembles that obtained at the one-loop order[l(|. Thus, higher order effects do not modify the 
qualitative peak-structure of the fermion spectrum at this coupling. There is, however, a quantitative difference 
for the central peak: While it is a (5-function shape at the origin in the one-loop calculation, which can be shown 
analytically, it has a finite width in the present result. We plot p+ at zero momentum in Fig. [3] from which we see 
that the width of the peak at the origin becomes broader and the peak height smaller as T increases. It is noted 
that the behavior of the peak height is consistent with that in the one-loop calculation where the residue of the 
^-function peak becomes smaller as T increases. The behavior of the peak height can be found from that of the 
self-energy at zero momentum shown in Fig. 3) We see that |Im£+(po> 0)| at the origin increases as T increases, which 
means that the decay amplitude of the fermion increases with T. This is contrasted to the imaginary part of the 
self-energy at one loop which vanishes at the origin independently of T*[To| . The reason of the vanishing at one loop 
is as follows. We can readily show that the fermion with the energy — m < po < m and the zero momentum decays 
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FIG. 5: The kinetic processes by the Landau damping contained in the imaginary part of the fermion self-energy at one loop. 
(The corresponding inverse processes are not shown.) The solid line represents the fermion and the dashed one the on-shell 
boson with the energy E v -k = — k) 2 + m 2 - The fermion with the energy k is also on-shell. 
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FIG. 6: The spectral functions p+(po,0)m for g 2 = 1 and T/m = 0.6, 1.2, 1.6,2 from top to bottom. 



through the Landau damping shown in Fig. [5] at the one-loop order. It is noted that the boson and the out-going 
fermion are on-shell in this figure. Then, in the left process of Fig. the boson momentum k must be infinity for 
Po = due to the energy-momentum conservation. Because this process and the inverse one give the statistical factor 
n(y/k 2 + to 2 ) + f(po + \/k 2 + to 2 ) in the amplitude, these processes are forbidden for the infinite k. Likewise, in the 
right process of Fig. [5j the momentum of the in-coming fermion must be infinity for po — and then the statistical 
factor vanishes as in the left process. Therefore, we have always ImS+(0, 0) — at one loop. Because ReE+(po,0) 
is an odd function by the charge conjugation invariance and thus vanishes at the origin at any order, the spectral 
function has a (5-function shape at the origin. In the SD approach, on the other hand, higher loop effects do not 
forbid the decay of the fermion with po = p — and thus Im£+(0, 0) can be finite, while ReE + (0, 0) = due to the 
property of the odd function as shown in the right panel of Fig. 2J This is the reason why there can exist a peak 
with a finite width at the origin in the SD approach. (In the right panel of Fig. [4] the line po/m — 1 is also plotted. 
Crossing points between this line and Re£+ give spectral peaks if ImE + /ReS+ is sufficiently small there.) We note, 
however, that it is nontrivial whether the three-peak structure appears in the SD approach because it would disappear 
if the value of Im£ + (0, 0) is quite large even for T/m < 2. The result shown in Fig. [2] indicates that the three-peak 
structure is stable against higher order corrections at this coupling. 

Next we examine the peaks in both sides which correspond to the normal quasifermion for pq > and the anti- 
plasmino for po < 0. In Fig. [51 we plot the spectral functions p+ for po > and p = obtained from the SDE and the 
one-loop calculation. Because we have a relation p + (po, 0) = p+(— po, 0), the anti-plasmino peak is also shown in this 
figure by replacing po with — po in the horizontal axis. At T/m = 0.6 and 1.2, the peak position and the width from 
the SDE are similar to those from the one-loop calculation. For higher T, the peak position and the width from the 
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FIG. 7: The fermion spectral functions p+m for T/m = 1 and g 2 = 0.4 (left), 2 (middle), 4 (right). The figure is clipped at 
p+m = 50. 

SDE are smaller and broader than those from the one loop, respectively. The broadening of the peak from the SDE 
could be understood from the fact that the probability of the boson emission and absorption from a fermion increases 
owing to multiple scatterings with bosons through the self-consistency condition The decrease of the peak position 
for higher T due to nonperturbative effects is a new result, for which we do not have intuitive explanation in the 
present. As seen in the next subsection, this behavior of the decreasing peak position is also seen when the coupling 
g becomes larger. 

B. Coupling dependence 
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FIG. 8: The imaginary part of the fermion self-energy E+/m for p — 0, T/m = 1. 

In this subsection, we study the coupling dependence of the fermion spectrum, focusing on the peak structure. In 
Fig. [7J we plot the fermion spectral function p+ for T/m — 1 and g 2 = 0.4,2,4. (For T/m = 1 and g 2 = 1, see 
Fig. [21) The fermion spectrum for g 2 = 0.4 is closer to the free spectrum than that for g 2 = 1 in Fig. [21 as it should 
be. For g 2 > 1, all the peaks at the low momentum region become broader and the height of the central peak falls 
rapidly. These properties are consistent with our previous work in which the fermion spectrum in a strong coupling 
gauge theory was investigated with the SDEQ: Because the SDE contains effects of multiple scatterings with bosons 
as mentioned in the previous subsection, the peak broadening for larger coupling could be understood from this point 
of view. The depression of the central peak is also understood from the multiple scattering effects: We plot the 
imaginary part of the fermion self-energy in Fig. [8] While the value at the origin is zero independently of the coupling 
at one loop as mentioned previously, its absolute value from the SDE gradually increases as the coupling increases. 
This means that higher order multiple scattering effects can bring about the decay of the fermion with po = p = 
and enhances as the coupling increases. 

Next we consider nonperturbative effects of the fermion spectrum by comparing with the result of the one-loop 
calculation as shown in Fig [51 At g 2 — 1, the peak positions and the widths are close to each other. At g 2 = 2,4, 
the widths are broader than those at one loop, which is due to the multiple scattering effects as mentioned above. 
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FIG. 9: The spectral functions p+(po, 0)m for T/m — 1 and g 2 = 1, 2, 4 from top to bottom. 




FIG. 10: The peak position of the normal quasifermion for p — and T/m = 1. 



The peak positions are smaller than the one-loop results. More specifically, we plot the coupling dependence of the 
peak position at po > in Fig. [101 We see that the peak position from the SDE is smaller than that from the one 
loop and clearly saturates around g 2 — 3. This behavior is consistent with that obtained in Ref.[7]. In the HTL 
approximation, the width of the corresponding peak vanishes and thus the peak position gives exactly the thermal 
mass which is proportional to T and the coupling. Our results suggest that a simple extrapolation of the results by 
the perturbation theory to the strong coupling region is not valid for the thermal mass. 



IV. CONCLUDING REMARKS AND OUTLOOK 



In this paper, we have investigated the fermion spectrum at finite temperature (T) in a chiral U(l)iX U(l)/j 
invariant linear sigma model employing coupled Schwinger-Dyson equations (SDE) for both fermion and boson with 
the ladder approximation. We formulate the SDE in the real time formalism and thus the spectrum in both the 
timelike and spacelike regions can be directly evaluated. To reduce the computation task, only the imaginary parts 
of the self-energies are evaluated with loop integrals, while the real parts are evaluated from the imaginary parts with 
the dispersion relation. 

We have assumed the massless fermion and the massive boson and focused on the T and coupling dependences of 
the fermion spectrum. The main reason is that it is pointed out from the one-loop evaluation in Yukawa models that 
the three-peak structure appears in the fermion spectrum at some T and coupling [ToL fllj . The essential points of the 
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existence of the three-peak structure found in Ref. [T(| are the massiveness of the boson and the Landau damping 
at finite T. One of the remaining issues is whether such a structure is stable against higher order corrections, in 
particular, for the strong coupling region. Our present result supports the existence of the three-peak structure even 
in the strong coupling region, although the strength of the central peak rapidly suppresses owing to multiple scattering 
effects with bosons as the coupling increases. The quantitative difference between our result and the result from the 
one-loop calculation is that while the central peak at the origin is of the (5-function shape at one loop, it has a finite 
width which becomes broader as the coupling and T increase in the Schwinger-Dyson approach. 

The present formalism can be readily applied to other models, such as gauge theories in the ladder approximation 
with the fixed gauge. Furthermore, because lattice QCD at high density is not realistic in the present, the extension 
to finite density systems is also important and, in fact, straightforward. 

Because the soft bosonic modes are expected to arise dynamically near but above the critical temperature of the 
second order transition, it would be also interesting to extend this formalism in order to study effects of the phase 
transition and temperature-dependent bosonic masses. This should be possible in future study within the presented 
dynamical framework. 
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